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Abstract 

We reexamine the scattering coefficients of shallow water waves blocked by a stationary counter 
current over an obstacle. By considering series of background flows, we show that the most relevant 
parameter is Unax, the maximal value of the ratio of the flow velocity over the speed of low frequency 
waves. For subcritical flows, i.e.. Unax < 1, there is no analogue Killing horizon and the mode 
amplification is strongly suppressed. Instead, when Unax 1.1, the amplification is enhanced at 
low frequency and the spectrum closely follows Hawking’s prediction. We further study subcritical 
flows close to that used in the Vancouver experiment. Our numerical analysis suggests that their 
observation of the “thermal nature of the mode conversion” is due to the relatively steep slope on 
the upstream side and the narrowness of the obstacle. 


1 Introduction 

In [1, 2], W. Unruh pointed out that one could observe the Hawking emission using analogue systems 
mimicking the scattering of light on a black hole metric. His argument rests on the possibility to 
observe the scattering of linear waves propagating against a transcritical counter flow. Since then several 
theoretical and experimental works have made this idea more concrete. Let us present here the most 
relevant results for the scattering of shallow water waves [3]. 

The first question that was addressed concerns the spectral properties of the analogue Hawking 
radiation taking short distance dispersion into account. It was found [2, 4, 5] that the spectrum is 
robust. That is, the relative deviations with respect to the Planck spectrum at the Hawking temperature 
k/2tt are linear in k/A, where n is the surface gravity of the hole, and A ^ k is the dispersive high 
momentum scale. In these works, the stationary background flow was assumed to be transcritical. This 
guarantees that long wave-length waves effectively propagate in an analogue space-time metric which 
possesses a Killing horizon. ^ 

The second question concerns the actual properties of the background flows that have been realized 
in experiments aiming at detecting the analogue Hawking effect. It appears that the background flows 
in water waves experiments [10, 11] (and probably also those using light in non-linear media [12, 13]) 
were not trans-critical. To understand the observations, one must thus start the analysis afresh and 
theoretically calculate the spectral properties in sub-critical flows. 

Following the procedure of [14], we study the scattering coefficients in a series of flows where Fmax, 
the maximal value of the ratio of the flow velocity over the speed of low frequency waves, ranges from 

^Subsequent works [6, 7, 8, 9] showed that the leading deviations scale as (k/A) x (for quartic dispersion), 

where D = KXiin/cn characterizes the spatial extension a;iin of the near horizon region where the gradients of the flow can 
be treated as constants, ch is the velocity of the low frequency waves evaluated on the horizon. Hence, to avoid large 
deviations one should work with D > 0.1, i.e., the flow must be clearly subcritical on one side and supercritical on the 
other side of the horizon. 
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1.25 (transcritical) down to 0.75 (subcritical). Our first aim is to display the key role played by i^max in 
dividing the behavior of scattering coefficients between trans- and sub-critical flows. For trans-critical 
flows with Fmax ^ l-l, we recover Hawking’s prediction up to small corrections. The scattering on 
subcritical flows instead does not follow simple laws. Keeping Fmax — 0.7 fixed, our second aim is to 
identify the subdominant role played by the length and the two slopes of the obstacle. For such flows, we 
show that changing any of these three parameters significantly affects the scattering coefficients. Finally, 
we compare our results to those obtained from the flow used in the Vancouver experiment [11]. Our 
findings suggest that the reported observation of the “thermal nature of the mode conversion” is related 
to the specific properties of their obstacle. 


2 Scattering coefficients in sub- and trans-critical flows 


We present the essential steps to compute the scattering coefficients of shallow water waves sent against 
a stationary flow. ^ More detailed explanations can be found in [14]. 

When the fluid is ideal, inviscid, and incompressible, the propagation of surface waves on an inhomo¬ 
geneous, laminar, two-dimensional flow is (approximatively) governed by [3, 14, 16, 17] 

[{dt + d^v{x)) {dt + v{x)d^) - igd^ tanh {-ih{x)da;)] x) = 0, (1) 


where the partial derivatives act on all factors on their right. Here v(x) is the horizontal flow velocity, 
h(x) the background fluid depth, and g the gravitational acceleration. The field x) is related to the 
linear variation of the free surface by Sh(t, x) = — ^ (dt + vdx) (j)(t^ x). Sh is the quantity measured in the 
experiments [10, 11, 18]. 

Since the background flow is stationary, we work at fixed frequency w, with stationary modes 
e~^‘^*(f>i^(x). As in [14], we work in a weak dispersive regime where Eq. (1) becomes 


{-iuj + dxv) (-ioj + vdx) - gdxhdx 


^dx (hdxf 


(j)ui{x) = 0. 


( 2 ) 


In the long wave length approximation, the quartic term is also neglected. Up to a conformal factor, 
the resulting equation has the form of a two dimensional d’Alembert equation in a curved metric given 
by = —(?dt^ + (dx — vdt)^, where c{x)‘^ = gh(x) is the speed of low-frequency waves. The metric 
possesses a Killing horizon if c(x) = If (a;)] at some point, i.e., if the Froude number F = \v/c\ crosses 1. 
Its surface gravity is given by k = c \dxF\ evaluated where F = \. 

In stationary flows, v(x) = J/h(x), where J is the conserved current. Hence the stream is completely 
characterized by h(x) and J. In what follows we study a series of flows parameterized by their minimum 
water depth dmin reached on top of the obstacle. When /imin < he = (J'^/g)^^^, Fmax = (hc/hxnin)^^^ is 
larger than 1 and the flow is transcritical. Such a flow possesses a black hole horizon on the upstream 
slope and a white hole horizon on the downstream one. Instead, when Umax < Ij the flow remains 
globally subcritical. Importantly, Eq. (2) offers a reliable approximation when the non-homogeneities 
of the flow are essentially localized in a region where F « 1. Hence Eq. (2) can be used to study the 
changes of the scattering coefficients when F^nax decreases from ~ 1.2 down to ^ 0.7. 

Scattering of incident counterpropagating waves 0^’*" was observed in [10, 11, 18]. For f > 0, these 
waves propagate to the left, as indicated by the exponentiated arrow. When scattered, they give rise to 
four outgoing waves: 


O-UJ (t>u 


,d,out 





( 3 ) 


Each outgoing wave is unambiguously identified by its wave-vector in an asymptotic region, see Fig. 1. 

^In this work, we only consider flows without undulation, that is, without the zero-frequency modulation of the free 
surface which is generally present in the downstream region. In Section III.B.3 of [14], we numerically found that the relative 
modifications of the scattering coefficients are of the order of Sh/h^in, where hmin is the minimum water height found on 
top of the obstacle, and 5h the amplitude of the undulation. For the flow of [11], we estimate that ~ O-l- Hence, 

the undulation should not play any significant role for this type of flows. The reader interested in the spectral modifications 
induced by undulations can consult Ref. [15]. Although that work deals with flowing atomic Bose condensates, we checked 
that its main results also apply to water waves. This was expected in virtue of the similarities of the Gross-Pitaevski and 
Korteweg-de Vries equations. 
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Figure 1: We show the dispersion relation associated with Eq. (1) in the fluid frame, VP' = gktscnhihk). 
The green continuous (purple dashed) curves describe the branches with positive (negative) values of Vi. 
The straight line gives fl = uj — vk evaluated in a subcritical flow to the right 0 < u < c = and for 
a fixed w > 0. The black dots indicate the four real roots k^. They are labelled by the coefficient of the 
corresponding outgoing wave of Eq. (3). Notice that only the last root lives on the negative branch. 


Starting from the left, the first root describes a dispersive short wave number mode propa¬ 

gating to the right. (By dispersive, we mean that this root exists only because the dispersion relation is 
not linear.) The second root describes a low wave number mode propagating to the left. When evaluated 
on the left (upstream) side, it describes the transmitted wave On the right side it describes the 

incident wave The third root describes the elastically scattered mode which also propagates to 

the right. As we shall see, this long wave length mode plays no significant role in the scattering. The 
fourth root instead is essential. It describes the other dispersive short wave number mode propagating 
to the right. Unlike the other modes, it carries a negative energy and has a negative norm.^ It has been 
complex conjugated so that all modes have a positive unit norm. 

Using these normalized modes, the scattering coefficients obey 

\a^\^ + \A^? + \A^\^-m^ = l. (5) 

Notice that one recovers the standard relation |q!^|^ — = 1 iff -I- ^ 1. In relativistic 

settings, when neglecting the gray body factor (the equivalent of |A„p), it is fulfilled and the radiation 
emitted from a Killing horizon follows a Planck law at the Hawking temperature: — 1)“^. 

Importantly, when Fj„ax > Ij this result is recovered from Eq. (2) in the dispersionless limit —>■ 0 
at fixed v{x) and c{x), where k is given by the expression below Eq. (2). As shown in [6], the leading 
deviations are proportional to h which here gives the dispersive short distance scale. 

2.1 Dependence on 

We now study the behavior of the four scattering coefficients when “removing” the Killing horizon by 
lowering Umax below 1. In Fig. 2 we show the series of 6 flows we shall use. These flows interpolate from 
Tjnax = 1-25 down to 0.75. The profiles of h{x) and F{x) = {hc/h{x))^P are given for a fixed value of the 

®The conserved scalar product of two solutions of Eq. (1) is given by [2] 

(4>i,(I> 2) = i J + vdx)(l >2 - + vd^)4il) dx. (4) 

Since is not positive definite, the scattering can lead to a mode amplification while preserving the norm. This 

mechanism is also known as over-reflection [19], and is at the root of the (analogue) Hawking effect. 
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current J = 0.126 m^/s. The intermediate flow with Fmax — 0.95 has been obtained by solving nonlinear 
hydrodynamical equations, following the method presented in App. A of [14]. The 5 other flows have 
been obtained from it by adding a constant value to the water depth. It should be noticed that their 
downstream slope (on the right) is steeper than their upstream one. 

h F 




Figure 2: We show the water height h{x) (left) and the corresponding Froude number F{x) = v{x)/c{x) 
(right) for 6 different flows to the right. The unit of the horizontal (and vertical for h) axis is the meter. 
The narrower profiles in dashed lines give our estimations of h and F for the flow used in the Vancouver 
experiment [11]. The horizontal dotted line on the right panel gives the critical value F = 1. We see that 
three upper flows are transcritical, whereas the other ones are subcritical. Notice that the undulation 
for the flow of [11] is not shown, see footnote 2. 

It should also be noticed that these flows radically differ from that of Ref. [11]. As can be seen in 
Fig. 2, the latter is much more narrow and has an higher slope in its upstream side, while it has typically 
the same slope on the downstream side. Its characteristic length scales are /imin ~ 8cm and /imax 20cm, 
and the corresponding current is J = 0.045m^/s. As we shall see in the second part of our analysis, its 
narrower character does affect the spectra, although less importantly than the value of Ffnax- ^ 

In Fig. 3, we represent the 4 coefficients of Eq. (3) obtained by numerically solving Eq. (2). On 
the horizontal axis, we use ln(a;/a;niax), where Wmax is the theoretical maximal possible frequency for a 
counterpropagating wave sent by the wave maker. For the present series, it ranges from 2.17 to 2.85Hz. ® 
On the vertical axis we show the logarithms of the norms of the four coefficients. From the two upper 
plots, giving respectively Inja^j] (left) and lnj/3^j (right), it is clear that the coefficients belong to two 
distinct classes depending on the trans- or sub-critical character of the flow. 

When the flow is sufficiently transcritical, the norms of aoj and grow as (Teff/w)^/^ for w —>■ 0. 
For the flow with the highest value of Fmax 1.25, the relative difference between Teff and the Hawking 
temperature k/2tt (calculated at the white hole horizon on the downstream side) is less than 13%. 
Further simulations show that Hawking’s predictions are recovered for a large range of frequencies when 
Emax ^ l-l- (For localized obstacles, the finite size of the transcritical region leads to a suppression of 
\Pui\ at ultra low frequencies not represented here, see [14].) These results are in agreement with those of 
Refs. [7, 9] which were obtained by considering transcritical flows with F{x) monotonically varying with 
X. The agreement can be understood on the basis that the emission spectrum is essentially fixed by the 
gradient on the downstream side. 

^Given the available information and the complexity of hydrodynamical flows, it is difficult to give a precise estimate 
of Enax for the experiment of Ref. [11]. The Bernouilli equation in the low-gradient approximation gives iTiax si 0.5. 
Taking into account the dispersive term from the Korteweg-de Vries equation gives iTiax ~ 0.65, close to the value used 
in [14]. Yet this value is rather sensitive to the parameterization of the obstacle one adopts. In brief, we estimate that 
0.6 < Enax ;S 0.75. The validity of this range is strongly supported by the good agreement with the observations of the 
transmission coefficient performed in [18]. In Fig. 2 and forthcoming simulations, we work with iTiax ~ 0.75. 

^Notice that Eq. (1) and the nonlinear equations of App. A of [14], are left invariant under the rescaling uj —> Acj, 
X —> \~^x, h —> \~^h, and J —> J, which leaves iTiax unchanged. This invariance is lost when taking into account 
capillary effects [20]. 
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Figure 3: The logarithms of the four scattering coefficients associated with the above series of flows are 
represented as functions of In . From the two upper plots one sees that the behavior of \auj\ (left) 

^max 

and 1/3^1 (right) principally depends on the trans- (the three upper curves) or the sub-critical character of 
the flow. In the left lower plot, for the three subcritical flows, one sees that the norm of the transmission 
coefficient approaches 1 (for decreasing values of w) near the critical frequency where the norm of a^) 
starts diminishing. The three vertical dashed lines indicate the values of Wmin, see Eq. (6). In the right 
lower plot, we see that the norm of the reflexion coefficient remains smaller than . 


Instead, for the three subcritical flows, \j3ui\ and \ai^\ vanish as for w —>■ 0 [14]. In addition, one 
notices a sudden decrease of ja^j for w smaller than a critical frequency, that we call Wmin- It is given 
by the double root of the dispersion relation evaluated on top of the obstacle. For near critical flows 
1 - Tmax < 1, one finds 

( 6 ) 

where Cmin = As can be seen in the left lower panel, the norm of the transmission coefficient 

\A^\ is close to 1 for a; < Wmim here indicated by three vertical lines. For uj < Wmin, the waves are 
thus essentially transmitted. This is easily explained by considering the characteristics of Eq. (2): the 
classical trajectories followed by counter propagating wave packets have no turning point for w < a;„iin- 
The large increase of transmission for oj < Wmin was clearly verified in a dedicated experiment [18]. On 
the right lower panel, for all flows, one observes that the norm of the “gray body coefficient” remains 
< Hence it plays no significant role in the scattering. 

In brief, we have shown that the scattering coefficients behave very differently in trans- and sub-sonic 
flows. For trans-critical flows with Emax ^1-1, the situation is clear as we recover Hawking’s prediction. 
As a consequence, since the temperature is hxed by the surface gravity n, the upstream slope and the 
size of the trans-critical region on top of the obstacle play no significant role. For the other flows and in 
particular for sub-critical ones, to our knowledge, the relevant parameters which determine the emission 
spectrum have not been identihed. 

2.2 Relevant parameters in snb-critical flows 

To guide our search, we present on the left panel of Fig. 4 the logarithms of the four scattering coefficients 
numerically evaluated for a flow similar to that used in the Vancouver experiment [11]. The characteristic 
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frequencies are Wmax 4.7Hz, and Wmin 1.8Hz. By comparison with the plots of Fig. 3, one first sees 
that the behavior of these coefficients belongs to the class of sub-critical flows. In particular, the decrease 
of \aui\ and the associated increase of for w < Wmin are both clearly visible. In addition, we also see 
that \a^\^ ^ l/3wP oc w for a; —?> 0, which is less visible from the two upper plots in Fig. 3 but which is 
common to all our numerical simulations when considering subcritical flows. 




Figure 4: On the left panel, we show the logarithm of the four scattering coefficients for the flow of the 
Vancouver experiment as functions of ln(a;/a;max)- The continuous, dashed, dashed-dotted, and dotted 
curves respectively correspond to |a^|, |/3^|, and \A^\, see Eq. (3). The decrease of la^^l and 

the corresponding increase of the transmission coefficient \Ai^\ are both clearly observed for uj < Wmin, 
indicated by a vertical dotted line. On the right panel, we represent the logarithm of the ratio R = 
\l3ijj/au)\^ for the Vancouver flow (in red dashed), for the three sub-critical flows, and that with highest 
Tmax of Fig. 3. The last curve is almost straight, as expected since |/3t.ap follows a Planck spectrum, 
whereas the curves for the three sub-critical flows show oscillations. 

Yet, we observe some interesting differences. First, the higher value of 1/3^^^ (with respect to that of 
the longer obstacle of Fig. 2 with the same value of Fmax and a similar slope on the downstream side) 
indicates that the higher slope of the Vancouver obstacle on its upstream side should play a significant 
role. Moreover, its narrower character should also explain why the decrease of the ja^jl is less pronounced 
than that of the longer obstacle. We verified these two conjectures by varying the slope of the upstream 
side and the length of the plateau on top of the obstacle, while keeping the value of F„iax and the 
downstream slope fixed, see Appendix A. In Appendix B we show profiles of perturbations of the free 
surface to illustrate the difficulties to extract the transmission coefficient at low frequencies. 

It is also interesting to study the behavior of the quantity that was observed in the experiment of [11], 
namely the logarithm of associated with the two dispersive waves, see Eq. (3). On the 

right panel of Eig. 4, as a function of uj (in Hz), we represent InA^; (numerically computed using Eq. (2)) 
for the Vancouver flow (in dashed), and for four flows of Fig. 3, namely the transcritical flow with highest 
Fmax and the three sub-critical flows. As noticed in [14], InA^, numerically computed for the Vancouver 
flow is almost linear, in agreement with what was observed in [11]. This is surprising since is highly 
suppressed, see the dashed curve on the left panel. In fact, Ini?^^ closely follows the one evaluated for the 
transcritical flow, which is there linear because \ j3uj\^ is Planckian, and — |/3wP = 1 is well satisfied. 
The important lesson is that the behavior of In R^^ is unable to differentiate between the radically distinct 
behaviors of and for these two flows. The same lesson applies to flows where F monotonically 
varies with x, see Fig. 6 of [14]. 
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Yet, it would be nice to understand why Inii^^ for the Vancouver flow is almost linear in w. ® On 
the right panel of Fig. 4, we notice that the sub-critical flow with the lowest value of Fmax — 0.75, 
comparable to that of the Vancouver flow, has a larger value of In when compared to the other two 
subcritical flows, even though its emission spectrum \j3u}\^ is smaller, see right upper panel in Fig. 3. 
This reveals that, when lowering the value of F„iax, decreases faster than so as to increase 

the slope of Ini?^;. We then notice that for very low frequency the slope of Ini?^; for the sub-critical 
flow agrees with that of the Vancouver flow before decreasing and showing oscillations. This suggests 
that a more linear behavior would be obtained if narrowing the obstacle, as the oscillations result from 
interferences between the scattering from its two slopes, which are here well separated. (In fact, such 
oscillations in the spectrum were already observed in [21] when considering flows with two scattering 
zones, see also the black hole laser for similar interfering effects [22].) In Appendix A, we verify that this 
is the case. In brief, we conjecture that the linear slope of Ini?^^ observed in [11] is due to the rather low 
value of Fmax ~ 0.75 and the relatively narrow character of the obstacle, which is small enough to avoid 
interferences between scattering on the upstream slope and the downstream one. 


3 Conclusions 

By studying the variations of the scattering coefficients along series of background flows, we made the 
following important observations. Firstly, we showed that the scattering coefficients principally depend 
on Fjnax) the maximal value of the Froude number reached on top of the obstacle. Secondly, Hawking’s 
prediction are only recovered for transcritical flows with i^max ^ 1-1. Thirdly, for subcritical flows, we 
showed that varying separately the length and each slope of the obstacle significantly affect the scattering 
coefficients. Fourthly, we observed that oscillations of Ini?,^ found for long obstacles are suppressed for 
narrower ones. As a final comment, we wish to point out that although derived in the context of water 
flows, these results should also apply (up to minor modifications) to other systems such as flowing atomic 
condensates [23], polariton systems [24], and non-linear optical media [25]. 
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Appendix A: Length of the obstacle and linearity of Ini? 

We study the behavior of the scattering coefficients in subcritical flows when widening the obstacle by 
increasing the length of the plateau on its top and by decreasing its upstream slope. In the upper left 
plot of Fig. 5, we show the profiles of F{x) of the three flows. We see that their value of F^ax — 0.72 
is unchanged, and that their downstream (right) slope is essentially the same. The narrowest obstacle 
(dark green) has a length and slopes close to those of the obstacle used in the Vancouver experiment. 


F 



Figure 5: For three flows with the same value of F)naxj and approximatively the same downstream (right) 
slope, we show the Froude number (top, left) as a function of x, the scattering coefficients In |aa;| (top, 
right) and ln|/3a,| (bottom, left) as functions of ln(a;/a;niax), and lni?tj (bottom, right) as a function of 
oj. The green curves are obtained with an obstacle similar to the one used in [11]. The brown curves 
correspond to an obstacle with a smaller ascending slope, and the red ones to a longer obstacle. From 
the lower plots we clearly see that reducing the length of the obstacle, and increasing the downstream 
slope, increase the power and the linearity of Ini?^, even though one gets oc w as w —>■ 0 for 
the three cases. 

In the upper right plot, we observe that the slope of ja^jj for lu < Wmin significantly increases when 
increasing the size of the obstacle. This is expected since the transition between wave-blocking (w > Wmin) 
and transmission is smoother for narrower obstacles. In the lower left plot, the 3 curves representing 
In |/3ij| establish that, although oc w as a; —>■ 0 for the three flows, the reduction of the length of the 
obstacle, and the increase of the upstream slope, both increase the mode mixing for small frequencies. 
From this we conclude that, for subcritical flows, the downstream slope (of the “white hole” horizon) 
alone cannot explain the emission spectrum. This is confirmed by examining the behavior of the quantity 
observed in [11], see the lower right plot where InF^j is shown as a function of uj. We notice that reducing 
the length of the upper part of the obstacle significantly increases the linearity of Ini? in w. This can be 
understood from the fact that the radius of curvature for the narrowest obstacle is of the same order as 
the dispersive length h^[^. We conjecture that the absence of a clear scale separation between the values 








of the inverse gradient of the flow and the dispersive length Zimin for the Vancouver obstacle explains the 
linearity of In . 

We have performed extra simulations (not represented) where we multiplied, and divided, by 2 each 
slope. We found that decreasing the upstream slope monotonically and significantly decreases the slope 
of Ini?^. Instead, decreasing the downstream slope slightly increases that of Ini?^,^, while increasing it 
produces oscillations in In , probably because interferences between the two sides are enhanced. These 
simulations confirm that the upstream slope plays a leading role for this type of obstacles. It would be 
important to see if these properties (obtained from Eq. (2)) are validated in forthcoming experiments. 


Appendix B: Profiles of perturbations of the free surface 

We illustrate the difficulty of extracting the transmission coefficient from measurements of the free 
surface at low frequency, which means very long wave lengths. In Fig. 6 we present profiles of Sh^iix), 
the linear variation of the free surface for the incoming counterpropagating mode of Eq. (3) evaluated 
at 4 different frequencies. To see the transmitted amplitude, Shuj{x) is evaluated at a time when its 
amplitude in the downstream region reaches its maximum value. We used the parameters of [II], for 


CO ^ 
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Figure 6: Profile of the free-surface variation for the incoming counter-propagating mode scattered 
on a flow similar to that used in Vancouver, and evaluated at four frequencies uj = 2a;niin (top, left), uj = 
1.2a;inin (top, right), w = O.Swmin (bottom, left), and ui = O.lwmin (bottom, right). The corresponding 
values of \A^\ are respectively 0.055, 0.36, 0.70, and 0.99. The mode is normalized by imposing that the 
amplitude of the transmitted mode is independent of w, which means that the incoming wave amplitude 
remains finite in the limit w —> 0. 
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which Wmin ~ 1-8 Hz. On the right incoming side, the profile is complicated as 4 waves interfere. On 
the left it is simple as there is only the transmitted wave. When uj is significantly larger than Wminj one 
easily sees that the transmitted mode has a very small amplitude. When decreasing uj close to or below 
Wniirn this amplitude grows and becomes comparable to that of the dispersive waves. However, as the 
corresponding wave-length becomes very large (it is already of a few meters for oj = Wmim and goes to 
infinity in the limit a; —)■ 0), it is not easily seen. In fact, on the lower right plot, we notice that the 
incoming mode raises the height of the free surface on the right side by about the same amount as the 
transmitted amplitude. 

On the lower right plot, we also see that the typical amplitude of the short wave length modulations 
of on the right side is of the same order as the transmitted amplitude on the left, as if reflection 
was still significant. However, for w = O.lwmin one has \A^\ ~ 0.99 while the coefficients \a^\ and \I3^\ 
in front of the dispersive waves are both much smaller than 1, as can be seen on the left panel of Fig. 4. 
The reason for this discrepancy is as follows. The kinematical relations between and \a^\ (|/3w|) 

on one side, and between and \A^\ on the other side, imply that 


^1/2 ’ 


(7) 


see the discussion after Eq. (30) in [14]. Hence, when the ratio goes to constant for 

a; —>■ 0, as predicted by our numerical simulations, \ai^\ and \j3uj\ both vanish as These features 

might explain why transmission below Wmin was not reported in [11]. 
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